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Abstract 

Analytic results for the conductance of a molecular wire attached to meso- 
scopic tubule leads are obtained. They permit to study linear transport in pres- 
ence of low dimensional leads in the whole range of parameters. In particular 
contact effects can be addressed in detail. By focusing on the specificity of the 
lead-wire contact, we show that the geometry of this hybrid system supports a 
mechanism of channel selection, which is a distinctive hallmark of the mesoscopic 
nature of the electrodes. 



The nanometer length domain, the fate of the transistor integration length, is un- 
avoidably rising a new variety of phenomena concerning the role of quantum effects on 
electronics at the molecular scale [l|]. Although molecular materials for electronics are 
already realized 0, [| f|, [|, real molecular scale electronic devices 0, 0, § still have 
to cope with challenges in utilization, synthesis, and assembly ||. The idea of molec- 
ular rectification proposed since 1974 |TD| found experimental evidence only 20 years 



later [11, 12]. Concerning theory, conventional methods employed for characterizing 



transport properties in microelectronic devices, such as the Boltzmann equation [|13 



can no longer be applied at the molecular scale. Here transport properties have to be 
calculated by using full quantum mechanical approaches. 

Electron transmission through molecular and supramolecular interfaces was already the 
object of intensive theoretical investigation in the last decade due to the electron trans- 
fer phenomena underlying the use of scanning tunneling microscopes. More recently, 
studies of transmission properties of a molecular junction contacted to metallic leads 
HI, |I5!| have intensified the interest in the basic mechanisms of conduction across molec- 



ular bridges. In a parallel development the use of carbon nanotube (CNT) networks 
has been the focus of intense experimental and theoretical activity as another promising 
direction for building blocks of molecular circuits ITB, [T7J . Carbon nanotubes are known 
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to exhibit a wealth of properties depending on their nano-scale diameter, orientation 
of graphene roll up, and whether they consist of a single cylindrical surface (single- 
wall) or many (multi-wall) [18|, [HJ. However if carbon nanotubes are attached to other 



materials, the characterization of contacts P0| , [2T| becomes a fundamental issue with 
respect to the employment of tubes as elements of molecular circuits. This could be the 
case when a carbon nanotube is attached to a single molecule or a molecular cluster 
with a privileged direction along the current flow, namely a molecular wire. 
The archetype of a molecular device can be viewed as a donor and acceptor lead coupled 
by a molecular wire bridge. In such systems the traditional picture of electron trans- 
fer between donor and acceptor species is re-read in terms of a novel view in which a 
molecule can bear an electric current [p^j . In the usual treatment of transport through 
molecular wires, the attached leads are represented by a continuum of free or quasi-free 
states, mimicking the presence of large reservoirs provided by bulky electrodes. How- 
ever such an assumption may become inadequate when considering leads with lateral 
dimensions of the order of the bridged molecule, such as carbon nanotubes. The lat- 
ter have been recently used for enhancing the resolution of scanning probe tips [^3 



In a similar configuration the fine structure of a twinned DNA molecule has been ob- 



served |24[] . One can argue that the presence of such mesoscopic leads strongly influences 
the conductance properties across the molecular bridge. 

This paper addresses the influence of the molecular wire-electrode contact on the con- 
ductance when the structure of electrodes plays an important role. In order to isolate 
contact effects, we treat analytically a minimal model where the tubular leads are de- 
scribed by a rolled square lattice and the molecule by a homogeneous linear chain. This 
system supports distinct transport properties depending on the number and strength 
of contacts between the molecular bridge and the interface as well as on the symme- 
try of the channel wave function transverse to the transport direction. The validity 
of these analytical findings are then tested for real tubular systems by treating CNT 
leads (figure [l].) via a numerical recursive Green's function approach. Such a study 
delivers additional insight to recent numerical results on carbon nanotubes/molecular 
wire hybrids |25 |. In particular, we demonstrate that the configuration when only one 



molecule-lead contact is activated gives rise to complex conductance spectra exhibiting 
quantum features of both the molecule and the electrodes; on the other hand multiple 
contacts provide a mechanism for transport channel selection leading to a scaling law 
for the conductance and allowing for its control. Channel selection also highlights the 
role of molecular resonant states by suppressing details assigned to the electrodes. 
The hamiltonian of the full system H = H tn bes + -f^wire + -^coupling reads 

H = E E w 

a=L,Ri w i rc n at n 'a 

- E («m L «r C +h.C.)- T ™R ( a m R < lrC + h.C.) 

m L <M c m R <A/ c 

where the matrix element „, = e" 5 n n i —J?L „, \ contains the on-site energy of each 
of the n wirc = 1, . . . , N chain-atoms, e wire , the orbital energy relative to that of the lead 
atoms, e L ' R , and 7 L ' R , <-y wirc ; and T are nearest neighbour hopping terms between atoms 
of the left or right leads, molecular bridge, and the bridge/lead interface, respectively. 
Note that nL,R is a two-dimensional coordinate spanning the tube lattice. Summations 
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Figure 1: Scheme of the molecular wire-carbon nanotube hybrid with single (bottom) 
and multiple (top) contacts. In this paper, on-site energies e Q=L > R > wirc are fixed to zero. 



over niL and run over interfacial end-atoms of the leads. In general, there are 
M such atomic positions, defining the perimeter of the tube ends, the number of 
hybridization contacts range between M c = 1 (SC) and M c = M (MC). The square 
lattice tubes (SLT) are obtained by imposing periodic boundary conditions on the 
longitudinal cuts parallel to the lattice bonds of length a. In the case of CNT, when 
the graphene honeycomb lattice is rolled along the lattice bonds, armchair single wall 
(£, £) nanotube are obtained. In this case M = 2£. 

In order to derive transport properties, we make use of the Landauer theory [EH] which 
relates the conductance of the system to an independent-electron scattering prob- 
lem |27|. The electron wavefunction is assumed to extend coherently across the device 



and the two-terminal, linear-response conductance at zero temperature, g, is simply 
proportional to the total transmittance for injected electrons T(E-p) at the Fermi energy 

9 = X T(Ef) - (2) 
The factor two accounts for spin degeneracy. The transmission function can be calcu- 
lated from the knowledge of the molecular energy levels, the nature and the geometry 
of the contacts. It is given by 



T(E) = l^ual 



Tr {SS*} 



(3) 
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where Jl,Jr are quantum numbers labeling open channels for transport which belong 
to mutually exclusive leads, in our case the two semi-infinite perfect nanotubes. The 
molecular system attached acts as a scatterer, and S is the corresponding quantum- 
mechanical scattering matrix. The quantity |<Sj U - R | 2 is the probability that a carrier 
coming from, say, left of the scatterer in the transversal mode j'l will be transmitted to 
the right in the transversal mode jr. The sum in (|3|) is restricted to transversal modes 
whose energy is smaller than Ep. 

In the linear regime one can write the zero-temperature transmission as 

T(E) = I^rJ 2 6 (E ~ E ih ) 5 (E - E lR ) , (4) 

«L,«R 

by considering the zeroth order states confined to the left and the right electrodes 
and | in), respectively. T is the corresponding transition operator whose form depends 
on the details of the confinement. Note that these confined states are not the scattering 
states labelled by j'l and jr in equation (||) which are the eigenstates of the exact 
hamiltonian. Equation (fD is basically derived within a transfer hamiltonian treatment 
a la Bardeen (weak coupling assumption) |[28]| . It has been initially developed by Caroli 
et al. [p9| , and more recently applied to molecular wires by Mujica et al. |30], |31|. The 



relationship between the scattering matrix and the transfer hamiltonian approaches 



has been extensively worked out in molecular systems [^, |32| showing de facto their 
equivalence. This fortunate fact let us to make use of the formal result of the Bardeen 
treatment in the broader applicability context of the Landauer approach (which holds 
also in the regime of perfect transmission). In fact the transmission function introduced 
in equation ([3|) is exactly the same as used in equation (Q). One can see this by writing 
down the Green's function matrix of the problem G _1 = G wirc 1 + El + Sr written 
in terms of the bare wire Green's function and the self-energy correction due to the 



presence of the leads. Making use of the Fisher-Lee relation |33| one can finally write 



T = Tr {A L (E)G(E)A R (E)&(E)} = 4A L n (E)A* NN (E) \G 1N (E)\ 2 , (5) 

where A a (E) = i(S a - S at ) (E+iO + ) contain only one non-zero element (A\ ± and 
for the left and right lead, respectively), due to the geometry. The rhs of equation (^) 
coincides with the formula provided within transfer hamiltonian schemes [30, 31]. The 
matrix element A L ( R > is the left (right) lead spectral density which is related to the 
semi-infinite lead Green's function matrix Gi ea d- It is minus the imaginary part of the 
lead self-energy (per spin) 

^=l,r = J2 r^GwKmi). (6) 

m a ,m' a 

Finally one has to calculate from the Green's function 

(G^" 1 )^ ={E + i0 + - H wirc ) nn , = (E + i0 + - er e ) <W - 7$£), (7) 

the Green's function matrix element Gin in equation (|5]). It is given as a matrix 
element referring to the two ^-atom-molecule ends and is, computationally, an N x TV 
matrix inversion. Since only the molecular-end on-site energies are perturbed by the 
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interaction with the leads via the self-energy S a , some general conclusions can be drawn 
without the need of an explicit computation of Gin- 
The Green's function Gin, in equation (§), reads 



/^wirc 



Gin — o- (8) 

(-[ y>L /^wircN /i yR riw'vce\ yL yR ( fi-w\re\ z 
I 1 ^11^11 H 1 ^NN^NN ) ^ll^NN V^IN I 



The interaction with the leads dresses via the self-energy E a the bare molecular wire 
Green's function element GJ§°. The latter can be calculated analytically in the case of 
an homogeneous wire (£™ ire = e wire , 7^ ire n / \ = 7 wirc ). In fact projecting on the N dimen- 
sional molecular wire basis, the determinant of the bare molecular Green's matrix (0) 
factorizes a dimensionless function of only the number of chain atoms, and of the ratio 
8 = (E — e wirc )/(27 wire ). This leads to a closed form for the molecular contribution 
in the conductance. Namely, one can easily check that Gj§ e = ^ wireAr 1 det (G wirc ) = 
7 wirc ~ 1 £(0)/£(iV), and GJl rc = G^ = 7 wirc " 1 £(iV - l)/£(iV), where the exact form of 
£ reads: 

/ y \N+1 , nV+1 

i{N) = (8 + VS^T) -[£- VS^IJ . 
After some algebra one finds that £ possesses the following recursive property, 

m = (9) 

which leads us to re-write equation (§) as 



— £2L_ = c(n) -(^- + ?m) e(iv - 1) + 3£L(jv - 2). 



This means that the inverse of the Green's function matrix element connecting left and 
right leads can be written as a sum of the terms, representing the inverse of the bare 
Green's function matrix elements for a wire of N, N — 1, and N — 2 atoms. In the 
limit of weak contact coupling the behavior of the Gin element is dominated by £ (N) 
leading to N transmission resonances in the conductance of unit height. Nevertheless, 
if the effective coupling between the molecule and the lead is much larger than 'y wire ; 
£(N — 2) will become the dominant term. As a consequence the conductance spectrum 
is effectively that of an (A^ — 2)-atomic wire 2"5fl . The resonant behavior inside the wire 



band (|£| < 1), and its modification due to the lead coupling is easily understood by 
writing the transmission in the following compact exact form valid for all A^ > 1 

T _ 4<5 2 sin 2 (0) 

" (sin(AT + 1)# - (5 2 - A 2 ) sm(N - - 2A sin N$) 2 + 4 (5 sin Ntf - X5 sin(N - 1)0) 

where the a = X — id = H/ 7 wire is the self energy of the leads (for simplicity assumed 
equal) normalized by the wire hopping. The parameter 0, defined by 



2 ' 



•& = cos" 1 8 = - In 8 VS " \ (11) 
2 8 + V8^T V ; 

189 



SLT, M = 10 (5,5) CNT, M = 10 




Figure 2: The normalized spectral density as a function of energy and active contacts 
is plotted for M = 10 possible atomic contacts available; on-site energies and hopping 
terms refer to a = L, R-leads. 



is real in the wire band giving rise to resonances for injected electrons matching the 
wire eigenenergies. Outside the wire band $ is pure imaginary (sin's are effectively 
sinh's) and the transmission as a power law dependence from energy and an exponential 
one from wire length, that is T ~ \2S\ for \£\ 3> 1 in agreement with previous 
results |34 . This analytic expression for the transmission provide the generalization of 
existing ones H3lj [32] , [35], |36| to nonvanishing real part of the self energies. The density 
of states M = — Im Tr {G} /n can also be written in a closed analytical form. One can 
in fact take advantage of the fact that, due to the wire homogeneity, all the diagonal 
elements but the first and the last coincide 



Gkk 



1 £(N - 1) - 2a£(N - 2) + a 2 £(N - 3) 
i< fc <7v ~ r** £(A0 - 2a£(iV - 1) + a^(N - 2) 



By using the parametrization ( |lT1) one can easily recast the DOS in the following com- 
pact form 

_ 1 iVsiniV^-2(iV-l)crsin(iV-l)^+(A^-2) ( 7 2 sin(iV-2)^ 
7T7 wirc m ~ sin(A^ + - 2a sin N$ + a 2 sin(N - 1)# _ ' 

The calculation of the spectral function S, or equivalently of the "surface" unperturbed 
lead Green's function is finally needed in getting the conductance. Due to the form of 
the interfacial coupling in our model the self-energy simplifies to 

^ = W S G lcad (m,m'), 

m,m' <M C 

where only surface terms enter in the sum over the states in the leads, and the effective 
coupling is defined as r eff = Yy/M c . 
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In Fig. the spectral function A = — ImE is plotted in the whole range of possible 
contacts between the SC and MC configuration. It is the result of the calculation of 
the surface Green function obtained analytically for SLT [07 , |38|| and numerically for 



CNT [J25|, leads. As a function of the number of contacts M c , the system interpolates 
between two different scenarios. In the MC case, it is effectively the spectral density 
of one-dimensional leads, obtained by Newns in his theory of chemisorption [^J . Only 
the channel without modulation in the transverse profile of its wavefunction contributes 
to transport |J7|, The two-dimensional character of the leads enters as an energy 
shift of 2 r y L,R for SLT and of r ) L,R for CNT, yielding an asymmetric density profile with 
respect to the atom on-site energy e L,R . In contrast, the SC spectrum is symmetric and 
richer due to the contribution of all available channels. Additional features characterize 



CNT leads [25, 3S 



In conclusion, we have shown that nanotube / molecular wire / nanotube hybrid systems 
exhibit novel features in the conductance. The conductance of a homogeneous molecular 
wire possesses an analytical form in the full regime of the wire parameters and allows 
for the insertions of a nonvanishing real self energy typical when considering nanotube 
leads. By tailoring the geometry and dimensionality of the contacts, it is possible to 
perform a channel selection. In the MC limit the conductance becomes independent of 
the topology of the tubular electrodes and transport is effectively one-dimensional. 
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